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ABSTRACT 



Context. Part of the very high energy y-ray radiation coming from extragalactic sources is absorbed through the pair production process on 
' ^ . the extragalactic background light photons. Extragalactic magnetic fields alter the trajectories of these cascade pairs and, in turn, convert cosmic 
' background photons to gamma-ray energies by inverse Compton scattering. These secondary photons can form an extended halo around bright 
' VHE sources. 

, Aims. We searched for an extended emission around the bright blazars Mrk 421 and Mrk 501 using the MAGIC telescope data. 
T-H 1 Methods. If extended emission is present, the angular distribution of reconstructed gamma-ray arrival directions around the source is broader than 
r* ' for a point-like source. In the analysis of a few tens of hours of observational data taken from Mrk 421 and Mrk 501 we used a newly developed 
' method that provides better angular resolution. This method is based on the usage of multidimensional decision trees. Comparing the measured 
^ _ shapes of angular distributions with those expected from a point-like source one can detect or constrain possible extended emission around the 
, source. We also studied the influence of different types of systematic errors on the shape of the distribution of reconstructed gamma-ray arrival 
directions for a point source. 

■ ■ ' Results. We present upper limits for an extended emission calculated for both sources for various source extensions and emission profiles. We 
obtain upper limits on the extended emission around the Mrk 421 (Mrk 501) on the level of < 5% (< 4% ) of the Crab Nebula flux above the 
energy threshold of 300 GeV. Using these results we discuss possible constraints on the extragalactic magnetic fields strength around a few times 
10-'^ G. 

Key words, gamma-ray astronomy - Cherenkov telescope - AGN halo 



1. Introduction 

Blazars are well-known extragalactic sources of Very High 
Energy (VHE) y-rays. Their radiation traverses over large dis- 
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tances through the extragalactic space filled with CMB (Cosmic 
Microwave Background) and EBL (Extragalactic Background 
Light) photons and could be absorbed via the pair production 
process. This effect can lead to a change in the spectral shape 
of observed rad iation and it can be used to con strain the EBL 
density (see e.g. IStecker et all (Il992l) . iMazin et al., (,2007i) ). 
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A possible extended emission aro und extragalactic source s 
of VHE y-rays was first discussed by lAharonian et al.l (Il994h . 
y-rays with energies larger than 10 TeV are strongly absorbed 
via a pair production process on EBL and CMB photons rel- 
atively close to the source. Secondary y-rays can be produced 
in a cascade initiated by those primary photons. Note that the 
redshift dependent energy for which the gamma-ray opacity is 

Ell 1 vary w ithin a factor of ~ 3 in di f ferent EB L models 
iiske et a ll 120041: IStecker et all 12006^ [Frances chini et al., 

2009t 



will be delayed and stretched to much longer time scales (up 
to ~ 10 ^ years at ~ 100 GeV energy for B ~ 1 0''- G , 
Neronov & Sem i koz (20071 |2009|) see also iPlagal ( Il995h : 



20081: .Gilmore et al 



I2OO8I: iPrimack et al 



Finke et al.i 



Magnetic fields are very non-uniform and their strength can 
vary within many orders of magnitude depending on the loca- 
tion in the large scale structure. So far they have been measured 
only i n the g alaxie s with strength ~ juG dKulsrud & Zweibel 
I2OO8I: iBeckl l2008h . the cores o f galaxy clusters (with in the 
inner 100 kpc) ~ several |/G dCarilli & Tavloit |2002|) . and 
near borders of few clusters with strength 10'^ - 10'^ G (on 
Mpc scales) ( IXu et al 1 12OO6I: iKronberg et al.L |2007|) . On larger 
distance scales the magnetic field strengths are not known 
but th ey have to be much weaker (see e.g. Ide Angelis et alJ 
(l2008h ). Moreover, the theoretical models predict very weak 
(B «; 10 G) extragalactic magnetic fields (EGMF) in voids 
in large scale structures, outside galax i es and galaxy clusters 
(Kronberg', 1994"; Grasso & Rubinsteinl l200lt IWidrowi I200I 
[Neronov & Semikoz, 2009). It is believed that voids make up 
a significant part of the space volume. 

Development of a secondary gamma-ray cascade depends 
on its original energy and the strength of surrounding magnetic 
field. Two cases can be considered. First, if the TeV source 
is located in an intense magnetic field region and the maxi- 
mal gamma-ray energy is large {E,„a_^ » 50 TeV), gamma-rays 
would produce e^e^ pairs directly near the source (within a few 
Mpc). These will be isotropized in the strong magnetic field 
whi ch might exist around t he source host galaxy or galaxy clus- 
ter (lAharonian et al.Lll994l) . 

Second, the mean free path for gamma-rays of moderate 
energies (E < 50 TeV) is longer, therefore they can tra- 
verse larger distances from the source. Therefore e^e^ pairs 
will be produ c ed in regions with much we a ker magneti c fields 
dPlagal 1 19951: iNeronov & Semikozl 120071: lElviv et all 120091: 
iDolag et all l2009h . In this case the deflections of the pair tra- 
jectories by EGMF are not large enough to make the secondary 
cascade emission isotropic. Instead the inverse Compton scat- 
tering of CMB photons by e^e' pairs will produce secondary 
y-rays with a slightly different direction with respect to the pri- 
mary photons. The re-direction of the cascade photons into the 
field of view of the telescope can lead to the appearance of an 
extended emissi on around the point-like source, even in the case 
B < 10 '^ G (INeronov & Semikozl l2007tlElvivet ail l2009t 
iDolag et"alll2009l) . In fact the extended emission would be pro- 
duced between the source and the observer. For a distant ob- 
server that extended emission will appear as superimposed onto 
the point source thus mimicking a halo. 

The secondary y-rays also can initiate cascades, provided 
that the optical depths in the EBL radiation field are still large 
enough. The energies of further generation of y-rays may be be- 
low the energy threshold of VHE y-ray instruments. 

Blazars are known t o be strongly var i able in particular at 
the y-r ay energies (e.g. lAharonian et al.l (l2007l) . lAlbert et alj 
(I2008al) 1. The extended emission component cannot follow 
the original time profile of the emission. This is because the 
secondary cascade photons do not propagate along the same 
path as the direct y-rays from the primary source. Instead it 



Dai et all (l2002l) : lMurase et alJ (12008*)). Thus the direct emission 



from the point source will be overlayed on the extended emis- 
sion component. The latter can constitute a part of the "quasi" 
constant, quiescent emission. 

The first attempt to detect extended emission around extra- 
galactic TeV y-ray sources was performed by the HEGRA in- 
strument for Mrk 501. It yielded only an upper limit of 5-10% 
of the Cra b Nebula flux (at energy > 1 TeV) on angular scales of 
0.5 to 1°, (lAharonian et al l l2001al) . 

In this paper we report on our search for extended emission 
of VHE y-rays from the bright blazars Mrk 421 and Mrk 501 us- 
ing the 17m diameter MAGIC Imaging Atmospheric Cherenkov 
Telescope (lACT). The details of the telescope, its performance 
and the standard analysis chain are described in lAlbert et al.l 
(I2008dl) . Due to the large size of the mirror dish and improved 
fight sensors, the MAGIC trigger threshold (~ 50 GeV) is ~ 2-3 
times lower compared to other operating lACTs. 

This paper is structured in the following way: In section 2. 
we describe the analysis method, which was used for search- 
ing for an extended emission. Also, we present a new method 
that improves the angular resolution. In section 3. we perform 
analyses of possible systematic effects which can mimic the ex- 
istence of an extended emission. Then, in section 4. we describe 
the data sample used for the analysis. In section 5. we present 
results of the analysis: the upper limits on the extended emission 
from Mrk 421 and Mrk 501. Finally, in section 6. we discuss 
how those upper limits constrain the strength of the EGMF. This 
is followed by a short discussion on future prospects and conclu- 
sions. 



2. Analysis method 

We parameterize t he images of a ir showers by using the so-called 
Hillas parameters ( lHillasl[T985h . The angular distance between 
the center of gravity of the image and the shower direction (so 
called DISP) is correlated with geometrical and timing proper- 
ties of the image. The arrival direction for every event can be 
estimated. The distribution of the squared angular distance be- 
tween the estimated and the true source position (0^) is narrow 
and has a peak at 6* = for a point-like y-ray source. In case of 
an extended source this distribution shall be broader. 

The DISP parameter is proportional to the ellipticity (1 - 
Width/Length) of the image. By including the dependence 
on the possible truncation of the image at the edge of the 
camera and also on the parameter Size (sum of total charge 
of an image), one can improve the precision of DISP (see 



or an image), one can irnprove tl 
iDomingo-Santamaria et al] ([2OOI ). 



2. 1 . The novel Random Forest DISP method 

For this study we developed a novel method for the DISP estima- 
tion that along with standard Hillas parameters (Width, Length, 
Size) includes the Leakage, the Time Gradient, and the depen- 
dence on the zenith angle of observations. Leakage is a measure 
of a truncation of images due to the camera edge effects. It is 
defined as the ratio of the charge in the last two rings of pixels 
in the camera to the total charge. The parameter Time Gradient 
is defined as the derivative of signal arrival time in pixels along 
the main axis of the image. For a given event the Time Gradient 
is strongly correlated with the impact parameter of the parent 
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Fig. 1. Comparison between the 0~ distributions of the y- 
ray excess obtained from 43h of Crab Nebula data with the 
parametrized DISP value (blue shaded area), and the RF DISP 
method (red triangles). 



shower Since for a given zenith angle and y-ray energy the 
DISP is a simple function of the impact parameter, one could 
expect to improve the DISP estimation, and as a result also the 
angular resolution by using the fast timing properties of the im- 
age. The MAGIC telescope's ultra-fast time response because 
of the parabolic reflector shape, special PMTs and 2 GSample/s 
FADC readout enabled us to include the Time Gradient parame- 
ter in our analysis thus enhancing the sensitivity of the telescope 
(lAliuetalll2009l) . 

In order to combine information from both the geometrical 
and timing properties of images in the most efficient way, we 
used multidimensional decision trees - the so-called Random 
Forest (RF) method. It is widel y used for the y/hadron separa- 
tion an d the energy estimation ( Bock et all. l2004t lAlbert et al.L 
120083) . A comparison between the novel RF method (hereafter 
RF DISP) and the standard parametrized DISP is presented in 
fig. [l] The RF DISP provides a substantially narrower dis- 
tribution and improves the angular resolution (defined as the 
40% containment radius for a point like source, equivalent to 
one standard deviation of a two dimensional gaussian distribu- 
tion) by ~ 20 - 30%. This improvement is due to the usage of a) 
Random forest method instead of simple parameterization and 
b) Time Gradient information. This enhances the telescope per- 
formance for the search for an extended emission. 

The shape and the width of a distribution for a point-like 
source depends on many factors, among them the energy of the 
showers. For a higher energy shower, due to the large number 
of particles in the shower maximum, one has a higher signal-to- 
noise ratio, and the resulting image has more precisely defined 
parameters. For this reason we have selected events with large 
S ize (the total measured charge of the image), which improves 
the precision of the reconstruction of the shower direction. In 
our analysis we used only showers with Size > 400 photoelec- 
trons which allows us to determine the 0^ with relatively high 
precision. This leads to an energy threshold (defined as the peak 
of the Monte Carlo simulated differential energy distribution) of 
300 GeV. 

Let us consider a situation where the excess observed from 
a hypothetical source is a mixture of a point source and a weak 
extended emission with a given profile. The cumulative (9^ distri- 
butions for point-like and extended sources are shown in fig. |2] 
Using this figure, the angular resolution for this analysis is esti- 
mated to be < o.r. 









0.02 


0.04 0.06 
theta^ [degree^] 


0.08 


0.1 




z \ 





















0.04 0.05 0.06 
theta^ [degree^] 



Fig. 2. Comparison of Monte Carlo cumulative 6^ distributions 
(upper panel) for a source with 80% point-like and 20% ex- 
tended emission and difference between them (lower panel). The 
characteristic radius of the extension is equal to 0.1° (red), 0.2° 
(green) or 0.3° (blue). The extended part of the emission is sim- 
ulated as a flat disc {dN/d&^ - const). The purely point-like 
source is shown as black line. A random mispointing up to 0.03° 
has been included in the simulations. 



2.2. Source extension 

If the characteristic extension 0e„ is much smaller than the tele- 
scope's PSF, the distributions of a purely point-like and a par- 
tially extended source are very similar. For larger extensions 
clear differences between the distribution shapes can be seen (see 

To investigate possible extended emission from blazars we 
adopted a method similar to the one used in lAharonian et al.l 
(I2001ah . We calculated the ratio of event rates in two ranges: 

N(0 <6<G ) 

f - -jjf^^g^gff, and compared them with the ratio calculated for 
a point-like source. The 0i and 02 values are calculated with the 
help of Monte Carlo simulations searching for the most signifi- 
cant difference between a purely point source and a source with 
an extended emission for a given size and profile of extension. 
For every considered flux, radius and profile of the extended 
emission we calculated the value of / and of the corresponding 
significance of the extension with MC simulations. 

To minimize the systematic errors in modeling a point source 
we have selected a data sample taken from the Crab Nebula 
(normalized to the Mrk 421 or Mrk 501 flux). The extension 
of the Crab Nebula in the VHE gamma-rays, as shown in 
lAharonian et a l. (2000) is below 0.025°, which makes it essen- 
tially a point-like source for MAGIC. Also the mean spectral 
slopes of both Mrk 501 and Mrk 421 in the considered data sam- 
ple are rather similar to that of the Crab Nebula. 



3. Analysis of the systematic effects 

Systematic effects can degrade the precision of the estimation of 
the arrival direction of y-rays. Some of those effects, if not taken 
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Fig. 3. Monte Carlo cumulative 0^ distributions for a point-like 
source for two different values of the optical PSF of the reflec- 
tor: a- = 0.036° (black), cr = 0.044° (red) (upper panel), and a 
difference between them (lower panel). 
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Fig. 4. Monte Carlo cumulative 9^ distributions for a point-like 
source with a different spectral index: -2.2 (blue), -2.4 (black) 
-2.6 (red) (upper panel), and a difference between them (lower 
panel). 



into account properly, can emulate an extended emission around 
the point-like source. Here we present a study of those effects. 

3.1. Optical PSF 

The optical point spread function (PSF) of the telescope smears 
out images. This increases Width and Length of any given im- 
age, thus changing its ellipticity. Therefore the DISP method will 
be affected and the 6^ distribution will become broader. Variation 
of the optical PSF will be reflected in the measured angular res- 
olution. The optical PSF varies across the field of view of the 
parabolic reflector of MAGIC. Also, varying gravitational loads 
of the mirror dish during observations produce small deviations 
of the optical PSF, which are corrected by active mirror control 
system of MAGIC. The largest observed variations of the PSF 
of MAGIC are in the range of 20 %. In fig. [3] we show the re- 
sult of a study of the influence of optical PSF on the shape of 
the 9^ distribution. An increase of the optical PSF by e.g. 20% is 
equivalent to < 2% admixture of an extended source with a char- 
acteristic extension radius of 0.2°. We conclude that this effect is 
negligible for our study. 

3.2. Spectral index 

As mentioned before, the 9^ estimation is more precise for higher 
energy y-rays. This means that a source with a harder spectrum 
will have a more peaked 9^ distribution than a source with a 
softer spectral index. 

In fig.|4]we present a comparison of cumulative 9^ distribu- 
tions for spectral indices -2.2, -2.4 and -2.6. The broadening of 
the 9~ distribution that corresponds to a steepening of the spectral 
index by 0.2 (a typical systematic error in spectral index deter- 
mination) is comparable to having a 3% admixture of extended 
emission with a characteristic extension of 0.2°. Compared to 
the statistical errors of this study, this effect is small. 




-0.04 -0.03 -0.02 -0.01 0.01 0.02 0.03 0.04 
Mispoint Az*sin(Zd) [degree] 

Fig. 5. Difference between the true and the estimated source po- 
sitions for Mrk 421 (red) and Crab Nebula (blue). Circles with 
radii 0.02° and 0.03° show the characteristic mispointing scale 



3.3. Mispointing 

A factor which limit the ability of lACTs to distinguish between 
point-like and extended sources is the possible mispointing of 
the telescope. Using a strong source like the Crab Nebula or 
Mrk 421 one can estimate the source position and then calcu- 
late the mispointing as a difference between the true and the 
estimated positions. As it is shown in fig. |5] the accuracy of 
the source position reconstruction is < 0.02 - 0.03° for these 
observations . This i s in agreement with our previous studies 
(lAlbert et al.L l2008dh . Those values are well below the angular 
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resolution of the telescope and the investigated extensions of the 
AGN halos. 

4. Observations 

MAGIC, consisting of two 17m diameter telescopes, is located 
on the Canary island of La Palma at the Roque de los Muchachos 
Observatory (at 2200 m.a.s.l). The data presented in this paper 
has been taken with the first MAGIC telescope. 

Mrk 421 and Mrk 501 are nearby blazars and well-known 
VHE y-ray sources. Their spectra have been measured up to 
~ 10 - 20 TeV (e.g. lAharonian et"al] (Il999bh . iKrennrich et"al] 
(l200l V At those energies absorption in EBL radiation fields be- 
comes important so one may expect an extended emission com- 
ponent due to an AGN halo. Both sources are strongly variable. 
They are being monitored by various VHE y-ray experiments for 
nearly 20 years. During this time both quiescent states with flux 
as low as ~ 0.15 C.U. (Crab unit) and giant flares with flux up to 
~ 10 C.U. were observed. 

During 1995-1999 (with the exception of short time flares) 
the flux registered by the HEGRA instrument from Mrk 421 was 
well below 1 C.U. In 2000 it increased to a lev el ~ 1 C.U., and 
it furth er increased up to ~ 2.5 C.U. in 2001 (lAharonian et al.l 
(l200l ). Observed flux could be described using power-law with 
a spectral index -2.4 and an exponential cut-off at the energy of 
3.4 TeV. In MAGIC observations of Mrk 42 1 performed between 
November 2004 and April 2005 the flux varied between 0.5 - 2 
C.U.. Also the cut-off en ergy in this time pe riod seems to be 
at lower value ~ 1 .4 TeV dAlbert et all (l2007ah ). In the low state 
the source seems to have a ste eper spectrum with a spectral i ndex 
~ -3 between 0.5 and 7 TeV (lAharonian et al.Ll2ml2003l) . 

Observations of Mrk 501 are equally interesting. Historically 
the strongest activity period for this source was in 1997. The 
mean flux observed by the HEGRA, the Whipple and the CAT 
telescopes in that year was 1.3-3 C.U. The same instruments 
observed Mrk 501 in 1998-1999 in a low state of ~ 0.15 C.U. 
The source became more active in 2000 resulting in a flux 
measured by the HEGRA and Whipple instruments at a level 
of 0.35 - 1.2 C.U . (see review of all those observations in 
lAlbert et al.1 (l2007bl) ). MAGIC observati ons of Mrk 501 in 200 5 
resulted in a mean flux of ~ 0.5 C.U. (lAlbert et al.l (l2007bl) ). 
while observations performed in 2006 show a low state a t ~ 0.2 
C.U. with a spectral index -2.8 (lAnderhub et all (l2009l) ). The 
spectra of Mrk 501 can be well d escribed by a (possibly curved) 
power-law (lAlbert et al.l (l2007bb ). 

For studing the possible extended VHE gamma-ray emission 
we selected recent MAGIC observational data from Mrk 421 and 
Mrk 501. To minimize systematic errors, only data taken at low 
(< 30°) zenith angles were used. 

Data from Mrk 501 have been taken in April/May 2008 in the 
so-called ON/OFF mode (where the source is in the center of the 
camera). After quality cuts 26h of ON data were selected. 50h 
of OFF data were used for the background estimation. Mrk 501 
was in a rather low state during the above-mentioned time period 
(mean flux in the sample ~ 15% C.U.). The entire spectrum can 
be well fitted with a single power law with a spectral index of 
-2.42 + 0.03stat ± 0.2syst. After correcting for t he abs orption due 
to the EBL, by using the iFranceschini et al] (HOOS") model, we 
obtain the source spectrum with an index of -2.24 + 0.03stat ± 

0.2syst. 

The Mrk 421 data was collected between December 2007 
- February 2009. After quality cuts 38h of data were selected. 
The data was taken in the so called wobble mode, where the 
source position was shifted 0.4° from the center of the camera. 
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Fig. 6. Comparison of the excess Q- distribution for 26h of data 
from Mrk 501 (red circles) and a point-like source (blue trian- 
gles, 17h of Crab Nebula data) (upper panel) and the diff'erence 
of both distributions (lower panel). 

The opposite (with respect to the camera center) position was 
used for the background estimation. Since the background data 
was taken simultaneously with the source data, the systematic 
errors are small. A disadvantage of this approach is the fact that 
if the source extension is as large as ~ 0.4°, the signal and the 
background regions start overlapping. This effect, which can re- 
duce the sensitivity, has been studied and included in the Monte 
Carlo simulations. In the analyzed data sample Mrk 421 was in a 
high state (~ 1.3 C.U.). Its spectrum in the wide energy range is 
best fitted with a flat (spectral index -2) power law with an expo- 
nential break point at the energy Ecut = 2.1 TeV). In the limitted 
energy range of interest (0.3 - 3 TeV) the spectrum from this 
data sample can be also fitted with an effective power law with a 
spectral index of -2.42 + 0.02stat ± 0.2syst. 

The 0^ distribution for a point like source was calculated by 
using the Crab Nebula data. We use 43h of wobble mode Crab 
Nebula data, taken between October 2007 and March 2009, and 
17h of ON mode data taken during December 2007 and January 
2008. This data-set can be described by a spectrum of -2.3 + 
0.02stat ± 0.2syst for the energy band (0.3 - 3 TeV). 

5. Results 

6*- distributions obtained for Mrk 501 and Mrk 421 are presented 
in fig.|6|and fig.H 

For both sources the 9~ distributions match with a corre- 
sponding point-like distribution. Using the first 12 bins, which 
contain most of the excess events we calculated ^^/«rfoy-7. 0/1 1 
(for Mrk 421) and 3.8/1 1 (for Mrk 501). 

The sensitivity for the detection of an extended emission de- 
pends on the extension size and profile. In our calculations we 
assumed a power-law profile of the emission (dN/dO oc (fi) with 
various steepness indices fi - 1,0,-1,-2. We performed the cal- 
culations for extended emission radii up to 0.5° for Mrk 501. In 
the case of Mrk 421, the data were taken in wobble mode which 
allowed the estimation of the background from the same data set. 
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Fig. 8. Upper limit on the flux {E > 300 GeV) of the extended emission from Mrk 501 (left figure) and Mrk 421 (right figure) in 
flie C.U. for diff'erent source profiles and extensions dN/dO <x < 9 < Ocm' fat disk (J3 - I) (blue triangles); 0.1° < < 6cut- 

/3 - -1 (red crosses). See more details in the text. 
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Fig. 7. Comparison of the excess 6^ distribution for 38h of data 
from Mrk 421 (red circles) and a point-like source (blue trian- 
gles, 43h of Crab Nebula data) (upper panel) and the difference 
of both distributions (lower panel). 



We calculated the upper limits for the extension radii up to 0.3°. 
In this case there is no overlap between the signal and the back- 
ground regions in the wobble mode observations. In fig. [8] we 
present upper limits on the flux of the extended emission calcu- 
lated for different extension radii and profiles for Mrk 501 (left 
figure) and Mrk 421 (right). 

Flux upper limits for an extended emission from Mrk 501 
for different energy thresholds are shown in fig.|9l Depending on 
the source profile the upper limit for the extension radius ~ 0.2° 
ranges from ~ 4%C.U. above 300 GeV to 6-7%C.U. above 1 
TeV. 




Fig. 9. Upper limit (in C.U.) on the flux of the extended emission 
from Mrk 501 for different energy thresholds, different source 
profiles and extensions dN/dO oc eP:0° < 6 < 0^^!- flat (J3 - 1) 
(dashed lines), 0.1° <0 < 0cut, P ^ -1 (solid lines). = 1000 
GeV (red), 600 GeV (green) and 300 GeV (blue). 



6. Discussion 

As can be seen in fig. [8] the most stringent upper limits for the 
halo search for Ey ^ 300 GeV are achieved for a source exten- 
sion of 0.2° - 0.25° for Mrk 501 (4% C.U.) and 0.2° - 0.3° for 
Mrk 421 (5% C.U.). We checked that the values of upper limits 
for the emission profile dNjde k 0^,0.1° < 6 < 0cut, ji - -2 ox 
f5 = Q are nearly the same as for the case off) - -1. 

For both sources, the best upper limits on the extended 
source flux are at a level of about 30% of the quiescent point 
source flux. Also, it is interesting to note that for an extension 
size > 0.2° one can observe only a marginal dependence on the 
emission profile. Extension sizes < 0.2°, becoming comparable 
to the telescope's PSF, provide worse upper limits. 

The non-detection of extended emission around Mrk 501 
and Mrk 421 in the 0.3-1 TeV energy range imposes restric- 
tions on the properties of the highest-energy y-ray emission from 
these sources and/or on the physical characteristics of the in- 
tergalactic medium around the sources. Cascade photons with 
an energy of Ey ^ 300 GeV result from absorption of the pri- 
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r i'/^ 

mary y-rays in the energy range Ey^^ ^ 20 LEy/0.3 TeV TeV 

which propagate over the distance Dy ^ 40/c |^£'y„/20 TeVj Mpc 
eronov &SemikozLl2009l) . Here k ~ 1 is a numerical factor 
which accounts for uncer tainty of the EBL models, wh i ch giv e 
0.6 < K < 2 5 forlKne iske et a ll (|2004|) : IStecker et alj ( l2006h : 
iPrimack et all ( l2009l) : lFranceschini d'ail ( l2008l) . 

Assuming that the EGMF streng th is much higher tha n 
B ~ 10 G, as in the halo model of lAharonian et alj d 19941) . 
the derived constraint on the extended source fluxes could be 
used to constrain the isotropic primary source power at energies 
above ^ 20 TeV. Such a constraint is especially interesting in 
the view of the recent discovery o f TeV y-ray emissi on from 

2006; Alber t et al.L 
Within the 



The non-detection of extended emission at 



300 GeV 



the nearby radio gala xies ( Aharonia n et al 



I2008bl: lAcciari et all 12009^ lAharonian et al 



12009) 



general AGN unification scheme ( Urrv et al. , 11991 ), the high en- 
ergy peaked BL Lacs (HBL), like Mrk 421 and Mrk 501, are 
believed to be the relativistically beamed versions of FR I type 
radio galaxies, like M87 and Cen A. Following the logic of the 
AGN unification scheme, the detection of TeV y-ray emission 
from M87 and Cen A indicates that HBLs produce both beamed 
and isotropic TeV emission. Of course, the existence of such 
an isotropic component of VHE y-ray emission from Mrk 421 
and/or Mrk 501 is difficult to verify, because the isotropic emis- 
sion would produce a much smaller contribution to the 0.1-1 TeV 
band point source emission. However, if the isotropic emission 
spectrum of Mrk 42 1 an d/or Mrk 501 extends, similarly to M87 
(lAharonian et all l2006l) . to energies Ey^ > 20 TeV, absorption 
of the isotropically emitted y-rays on the EBL leads to the pro- 
duction of an extended emission halo a round HBL. The detec- 
tion of an extended halo as discussed bv lAharonian et al ] (11994 
would, therefore, provide a direct evidence for the existence of 
an isotropic multi-TeV emission from HBLs. 

In order to derive constraints on the isotropic emission from 
Mrk 501 and Mrk 421 from the limits on the extended emission 
flux, one needs to estimate the fraction of the halo flux within 
the measurement region of the radius 0.1° -0.5°. The observable 
angular size of the halos around Mrk 421 and Mrk 501 (both at 
the redshifts z - 0.03 and distance D ~ 150 Mpc) at the energy 
Ey ^ 300 GeV is expected to be ~ Dy/D ^ 1 5°/c. Assuming a 
surfac e brightness profile dNy/dQ. ~ 1/^? like in lAharonian et al.l 
(fT99l . one could find that the region 6 < 0.5° around the source 
contains ~ 3% of the halo emission. The halo is expected to be 
more compact at 1 TeV, ^ 8.3°/c, so that the central 6 < 0.5° 
region contains ~ 6% of the halo flux. This means that the total 
isotropic luminosity of Mrk 421 and Mrk 501 is limited to be 
less than LhaioC^yo > 20 TeV) - 1.3/cLcrab(£' > 0.3 TeV) and 
i-haioC^ro > 36 TeV) ^ 1.3/fLci-ab(£' > 1 TeV). It is clear that this 
limit is not very restrictive. Our analysis shows that the isotropic 
luminosity of HBLs could be more efficiently constrained via the 
search for an extended emission from a more distant source (for 
which the angular size of extended halos is smaller) and/or via 
observations with wider field of view instruments. 

If the EGMF is much weaker than 10 G as already men- 
tioned, deflections of cascade e^e~ pairs are not strong enough to 
isotropize cascade y-ray emission. In this case no isotropically 
emitting halo around the point source is formed. The extended 
cascade source should appear more compact depending on the 
EGMF strength. Non-detection of extended emission constrains 
possible range of EGMF strength. In this short discussion we 
mention only limits on the EGMF with large correlation length. 
These limits could be extended to the case of arbitrary corre- 
lation length of the EGMF in a straightf orward way, using the 
formalism o f|N eronov &Semikozl(l2009h . 



could impose a bound on the EGMF only in the case when 
the time-averaged primary (beamed) source emission spectrum 
extends to the energies above Ey„ > 20 TeV (see above). 
Unfortunately, the specific of observations in the TeV band and 
extreme variability of TeV-emitting blazars do not allow one to 
derive the time averaged spectra of the sources. 

In the case of Mrk 421, a high energy cut-off in the spec- 
i Tum at Ecut 1 ~ 5 TeV has been repeatedly reporte d 
dKrennrich et al.Ll2001h (see however. iKonopelko etal] (l2008b ). 
If such a cut-off is intrinsic and present in the time-averaged 
spectrum, the source luminosity at 20 TeV is expected to be a 
factor of > 15 lower than the luminosity at 300 GeV (assuming 
that the intrinsic power law photon index is F ^ -1.7, close to 
the one measured by Fermi dAbdo et aUl2009l) ). In this case the 
flux of the extended cascade emission at the energies > 300 GeV 
is a factor > 15 lower than the point source flux. This is consis- 
tent with the upper limits derived above. In the scenario of in- 
trinsic cut-off in the Mrk 421 energy spectrum the upper bound 
on extended emission around Mrk 421, derived from MAGIC 
observations, does not constrain the strength of the EGMF. 

On the contrary, the time-averaged spectrum of Mrk 501 ex- 
tends, most probably to much higher-energies. No intrinsic high 
energy cut-offs in the low or high activity state of the source were 
reported. This implies that the source luminosity at the energies 
above 20 TeV can be of the same order as the luminosity at 300 
GeV. At the same time, the bound on the extended source flux at 
a level of ~ 0.04 C.U. at 300 GeV is by a factor ^ 4 lower than 
the point source flux at the same energy. If the intrinsic source 
flux at ~20 TeV energy will prove to be higher than one forth 
of the flux at 300 GeV, the bound on the extended source flux, 
derived from MAGIC observations, might impose constraints on 
the strength of the EGMF within the region D ~ Dy ~ AQk Mpc 
around Mrk 501. In fact, the spectrum of Mrk 501 from the an- 
aly zed data sample af t er cor rection for the absorption by using 
the lFranceschini et al.l (l2008l) model (which provides a relatively 
low level of absorption) has a rather hard spectral index of -2.24. 
This yields only a factor (20 TeV/0.3 TeV)"-^'* = 2.7 decrease in 
the SED from the energy of 0.3 to 20 TeV. However both the un- 
certainties of the EBL absorption models and the mere fact that 
Mrk 501 is known to be variable tells us that the time-averaged 
spectrum over a longer periods of time may be found, for exam- 
ple, softer than what has been observed for the used data sam- 
ple. Measurement of the intrinsic time-averaged source flux at 
20 TeV will be possible only after precise measurements of the 
EBL in the mid-infrared and regular monitoring of the source on 
year(s) time scale. Therefore we give below only a qualitative 
estimate of the range of magnetic fields which might be con- 
strained by the Mrk 501 data. 

Assuming that the correlation length of the EGMF is 
much larger than the inverse Compton energy loss distance 



of electrons with energies E^ 



EyJ2 



10 TeV, D, 



30 10 TeV] ' kpc, one can find that the size of the ex- 
tended source around Mrk 501 is expected to be ©ext - 

0.4° [r/3.5]"' [£y/300 GeV]"' [fi/lO-''* g] where t = D/Dy is 
the optical depth for the primary y-rays with respect to the pair 
production on the EBL dNeronov & Semikol. l2009h . The non- 
detection of an extended source with 0.15° < ©ext < 0.5° at an 
energy 300 GeV might constrain EGMF with a strength in the 
range 4 x 10"'^ < B < 1.3 x 10"'"* G. A significant secondary 
cascade emission at 1 TeV will be possible if there is no cut-off in 
the primary spectrum below 40 TeV. In this case non-detection of 
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the extened emission at 1 TeV can be used for excluding EGMF 
strengths up to B 8 x 10 G 

It is very desirable to increase the sensitivity of the halo 
search. An efficient way to enhance the sensitivity will be to use 
much longer observations, albeit if the systematics (e.g. the spot 
size, pointing accuracy) is well under control. For example if 
one would use 400h of on-source time (during 1 year, accept- 
able zenith angles) then a simple scaling as a square root of time 
will provide a sensitivity on the level of 1% C.U. 

Better sensitivity for a halo search could be achieved when 
the observed source is in a low emission state. In the high emis- 
sion state the extended component remains constant, while the 
point-like emission provides additional background, hence de- 
teriorating the sensitivity. This effect is especially important for 
the sources that have an extension comparable to the angular res- 
olution of the telescope. 

Improving the angular resolution of the instrument is an ef- 
ficient way of enhance the sensitivity for the halo search, es- 
pecially when the emission radius is small. It would help to 
disentangle the extended component from the primary, point- 
like component down to lower extensions. For example with the 
planned CTA-like telescope array one could ach ieve an angular 
resolution in the range of 2 arcmin (lHillasllI989l) . Accordingly 3 
times smaller halo extension compared to the current study could 
be probed. 

By measuring with an instrument one order of magnitude 
more sensitive than MAGIC (like CTA) it seems realistic to 
achieve a sensitivity of ~ 0. 1 % Crab for the halo search. This is 
true for sources with an extension > 0.1 - 0.2°. An observation 
time of a few hundred hours will be necessary for such a study. 

7. Conclusions 

We have developed a novel method based on the Random Forest 
algorithm for the estimation of the source position. This method 
improves the angular resolution by ~ 20 - 30%. 

In the used data samples consisting of 26h observations of 
Mrk 501 and 38h of Mrk 421 no extended emission has been 
detected around these sources. Our study showed that if there 
is an extended emission around Mrk 501, then its flux is < 4% 
C.U. (see fig. [8] left), for the analysis energy threshold of 300 
Ge V. For Mrk 42 1 the upper limits are less stringent, because the 
source was in a high emission state, which causes the tail of the 
point-like emission to extend into the halo region thus creating 
additional background. The constraint on the extended emission 
flux from Mrk 421 is < 5% C.U. (see fig.lright). 

We analyzed different types of systematic errors connected 
with the observations of extended sources. We found them to be 
negligible compared to the upper limits obtained on the flux of 
the extended emission component of Mrk 42 1 and Mrk 501. 

With the second telescope, the MAGIC system can be oper- 
ated in the stereo mode. By combining simultaneous information 
from both telescopes, the angular resolution will be improved. 
This will allow us to perform a search for the AGN halo with a 
better sensitivity. The stereo observations improve the sensitivity 
2.5 times at lower energies (~ 100 GeV). Moreover, the better 
signal-to-noise ratio reduces the possible bias due to systematic 
errors. It will allow us to perform a sensitive search for halos at 
lower energies, thus extending the range of the EGMF strengths 
investigated with this method. 

A future y-ray telescope project - the Cherenkov Telescope 
Array is aiming for a much better angular resolution. In this case 
lower halo extensions can be investigated with an improved sen- 
sitivity. 
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